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Abstract 

The model consists of a signal process X which is a general Brownian diffusion process and 
an observation process Y, also a diffusion process, which is supposed to be correlated to the 
signal process. We suppose that the process Y is observed from time to s > at discrete times 
and aim to estimate, conditionally on these observations, the probability that the non-observed 
process X crosses a fixed barrier after a given time t > s. We formulate this problem as a usual 
nonlinear filtering problem and use optimal quantization and Monte Carlo simulations techniques 
to estimate the involved quantities. 

1 Introduction 

We consider in this work a nonlinear filtering model where the signal process X and the observation 
process Y evolve following the stochastic differential equations: 

(dX t = b(X t ,t)dt + a{X t ,t)dW t , X = x , 

\dY t = h(Y t , X t , t)dt + v(Y t , t)dW t + 6(Y U t)dW t , Y = y . 

In these equations, W and W are two independent standard real valued Brownian motions. We 
suppose that the functions b, a, h, v, and 5 are Lipschitz and that, for every (x, t) e]0, +oo) 2 , 
6(x, t) > 0, v(x, t) > and a(x, t) > 0. 
Let a be a real number such that < a < xq and let 

7f = mf{u >0,X U < a} 

be the first hitting time of the barrier a by the signal process. As usually, we consider that inf = +oo. 
Our aim is to estimate the distribution of the conditional hitting time 

P(rf>t|jf) (2) 
for t > s > and where (J r Y)t>o is the filtration generated by the observation process Y: 

Tj = a(Y u ,u < s). 

More generally, we shall denote by (Tf)t>o the filtration generated by the process Z. 

Such a problem arises for example in credit risk when modeling a credit event in a structural 
model as the first hitting time of a barrier a by the firm value process X. Investors are supposed to 
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have no access to the true value of the firm but only to the observation process Y, which is correlated 
to the value of the firm (see e.g. [3j |U). We will typically suppose that we observe the process Y 
at regular discrete times to = < t\ < . . . < t m = s over the time interval [0, s] and intend for 
estimating the quantity 

P(r a x >t\Y t0 ,...,Y tm ) 

for every t > s. Note that, if t < s, then, applying the Markov property to the diffusion Y, the 
computations boil down to the case s = t. In [2], the quantity 



P inf X u>a \Y t0 ,...,Y tri 

\u€[s,t] 

has been estimated by a hybrid Monte Carlo-Optimal quantization method in the case where the 
observation process dynamics is given by: 

dY t = Y t (h(X t )dt + u(t)dW t + 5{t)dW t ) , Y = y , 

where v and S are deterministic functions. However, the approach used in the previous work does not 
apply to our framework because we want to compute the conditional distribution of a function of the 
whole trajectory of the signal process from to t given the observations from to s, with s < t. 



Example 1. A particular case of model £[]) one may consider is the following "Black-Scholes" case: 

(3) 



dX t = X t (ndt + ad W t ) , X = x , 

dY t = Y t {rdt + vdW t + SdWt), Y = y , 



so that 



or 



dY t v dXt , v s , . ~ 

—± = -—± + (r- f i-)dt + 5dW t , (4) 
Y t a X t a ! 

Observe that setting r = fi and a = v yields 

dY t dX t ~ 
- = — + 5dW t , 

meaning that the return on Y is the return on X affected by a noise (see e.g. Q). 

Of course, in this case, we may compute theoretically the expression (0 by noticing that: 

t* = inf{u >0,X U < a} 

= inf | it > 0, 5W U - 5 —u > ln(y n ) + In - ln(a) 

and that, conditionally to Fj , the process (w u , u < s^j has the same law as 

f B u + ~~o~~nr? ( ln (^«) + in ( — V ^4^ u ) ' u < s 



where B is a standard Brownian motion independent from the process Y. This follows from the fact 
that, since W and W are two independent Brownian motion, so are 

„ aW + SW A ~ SW-aW 
B = — ; and B = — , 

yV 2 + 5 2 V^TP 
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and from the relations 



1 A /^x , f x o\ a 2 + 5 2 \ J ~ SBt-aBt 



= (ln(Y t ) + In ( — ) + t ) 



Therefore, in this particular setting, the problem boils down to the computation of the first passage 
time of a Brownian motion to a curved boundary. We refer to Q for some similar, and far more 
general, considerations. 

Example 2. Another particularly simple case is given when both the signal and the observation pro- 
cess evolve following the Ornstein-Uhlenbeck dynamics: 



dX t = \(6 - X t )dt + adW t , X = x , 

dY t = \(9 -Y t )dt + adW t + SdW t , Y = y . 



(6) 



In this case, we have : 



so that 



X t =Y t - (y - x Q )e~ xt - 5e~ xt [ e Xu dW u , 

Jo 

t£ = inf |u > 0, Y u - a - (y - x ) e - Xu < 5e~ Xu e Xv dW v 



and, conditionally to , the process ye Xu j e Xv dW v , u < sj has the same law as 

° - Se Xu _ Xu Qn , u , 



u u + -^rnv [ Y u - - o{i - e~* u ) ),u< s 



where U is an Ornstein-Uhlenbeck process with parameters A and started from 0, and independent 
from Y. This follows from Knight's representation theorem combined with the same ideas as above. 

The rest of the paper is organized as follows: in Section[2j we state and prove the main theorems, 
i.e. we give an approximation of the expectation Q when X is replaced by its continuous Euler 
scheme X, see especially Subsection 12.31 Then, in Section 3, we introduce some numerical tools in 
order to compute the quantities involved. We finally conclude the paper by a few simulations. 

2 Estimation of the conditional survival probability 

2.1 Preliminaries results 

To deal with the computation of the conditional hitting time (0, we introduce the first hitting time 
from time s: 

t {s) =inf{u > s,X u < a}. 
Define furthermore the filtration (Gt)t>o by, for every t > 0, 

Gt = <r{(W u ,W u ),u<t}. (7) 

We have the following result. 



3 



Lemma 2.1. For every t > s, 
with 



E 



l {T x >s} F(s } t,X s )\T t 



Y 



(8) 



Proof. If t > s, one has 
so that 



P(r* > t\ J 7 , 



F{s,t,X s ) = P ( inf X u > a|X 

v wG[s,t] 



{^>t} = {r^>s}n{r^>t}, 

= P({r* >-}n{rW>t}|jT) 



E 



E 



l{T^> S } IL {r(»)>i}l^ 



where the a-algebra Q s is defined in ©. Hence, using the fact that W is independent from W and 
the Markov property of the process X: 



1 



{t*>s} 



E 



1 



{r( s )>t}l^s 



Y 



1 {ri c >s} E 



1\ I T w 



l {T x >s} P(r( s )>t|J-f)|J- s y 



l {T x >a} P( inf X u >a|X a )|j; 



This shows assertion ©. 



□ 



Note that in general, there is no closed-form expression for computing P ( inf X u > a\X s ) , except 

u£[s,t] 

in some very special cases, such as Brownian motion or Bessel processes... We shall therefore need 
an approximation of this expression, which is the purpose of the next subsection. 



2.2 The continuous Euler scheme 

Let us denote by X the continuous Euler scheme process associated to the signal process X and 
defined by 

X s = X^ + b(Xs,s){s -s) + a(Xs,s)(W s - W^, X = x , 

with s = tk if s G [tfc, tfc + i). To estimate the distribution of the conditional hitting time given in (O, 
consider that we observe the process Y at discrete and regularly spaced times: to,t\, . . . , t m , with 
= to < ■ ■ ■ < t m = s < t m+ i < • • • < t n = t. In our general model (see Equation CO), the 
discrete time observation processes {X tk , k = 0, . . . , n} and {Y tk , k = 0, . . . , m} are obtained from 
Euler scheme as: 

*t k+1 = X tk + b(X tk ,t k )A k + a(X tk ,t k )(W tk+1 - W tk ) 

Yt k+1 = Y tk + h(Y tk ,X tk ,t k )A k + u(Y tk) t k )(W tk+1 - W tk ) + 5(Y tk ,t k )(W tk+1 - W tk ) 

where k = 0, . . . , m — 1 for the observation process, k going till n — 1 for the signal process and 
A k = t k+ \ — t k . Note that if t > s, the number of discretization steps over [0, s] may differ from the 
number of discretization steps over [s, t] so that we choose 

t k = — for k = 0, . . . , m 

t k = s + k( ^~^ for k = m + 1, . . . , n. 
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Supposing that we have observed the trajectory (Y to , . . . , Y tm ) of the observation process Y we aim 
to estimate 

P(r* >t n \Y t0 ,...,Y tm ) 

by 

P(r* >t n \Y t0 ,...,Y tm ), (10) 

where 

= inf{u > 0,X U < a}. 

To this end, we will use a useful result called the regular Brownian bridge method. This result recalled 
below allows us to compute the distribution of the minimum (or the maximum) of the continuous 
Euler scheme X of the process X over the time interval [0,t], given its values at discrete and regular 
time observation points = to < t\ < ■ ■ ■ < t n = t (see, e.g. J6JI3). 

Lemma 2.2. Let {Xt)t>o be a diffusion process with dynamics given by 

X t = x+ / b(X u ,u)du + / a(X u ,u)dW u 
Jo Jo 

and let (Xt)t>o be its associated continuous Euler process. Then, the following equality in law holds: 
£{ u ^ t] Xu\X tk =x k ,k = 0,... ,») =£(^^ i (^3+^ fc) )- 1 (A*)) > (11) 



where (A k ) k= o ... n ~\ are i.i.d. random variables uniformly distributed over the unit interval and 

(' 

by 



jjXk,Xk+i \ is the inverse function of the conditional cumulative function H^' X 2? 1 v defined 



H *k,xk+i (u) . = f exp ( - Aka 2 2 {xkJk) {u - x k )(u - x k+1 y) ifu < mm(x k , x k+1 ) ^ 
A k a {x k ,t k ) | j otherwise. 

In the following, we shall replace the expression a(x k , t k ) in the expression of H by a k and rather 
write: H x A k ' x % +1 . 

Short proof. Observe first that, conditionally to {X t , = x k ; k = 0, • • • , n}, the random variables 



min X u ; k = 0, . . . , n — 1 > are mutually independent, thanks to the independent increments 

ue[t k ,t k+1 ] J 

property of Brownian motion. Then, it suffices to notice that computing the law of min X u 

«e[tfc,tk+i] 

conditionally to {X tk = x k ,X tk+1 = x k+ \} amounts to computing the law of the hitting times of the 
bridge of a Brownian motion with drift. This law is well-known to be independent from the drift (i.e. 
from the function b here) and to be given by an expression such as (fT2l) . □ 

In the rest of the paper, the function 1 — fj^' x ^ +1 ( a ) w jn be often used. We shall denote it by 



gXfc.asfc+1/ x that 



nx k ,x k+l , s A / 2(x k - a)(x k+1 - a)\\ 

fcf7 fc (a) = I 1 " 6XP { A k aHx k ,t k ) ) ) ^>-- (13) 

Now, in our set-up, we shall apply the Brownian bridge method to obtain the following lemma, which 
is taken from O : 
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Lemma 2.3. We have 



P( inf X u >a\(X tk ) k =0,...,n 

\UG [tm )tn] 



n-l 

n 

k=m 



<'„f" + '(a). 



(14) 



Proof. It follows from Lemma [2721 that 



P( inf X u > a\(X tk ) k =o,...,m 
\ue[t m ,t„] 



min 1 H A k «l 



A fe cr fe 



(A*) > a . 



Since the function 2 +1 (0 is non-decreasing and the A^'s are i.i.d. uniformly distributed random 
variables, this reduces to: 



min X u > a|(A\) fc=0 ,..., m = J] P f A* > H ^ 



U^z [tm ,tn 



n-l 



k=m 
n-l 



T Xt k' Xt k + l , 



k=m 



n 



(a) 



This completes the proof. 



□ 



2.3 Main theorem 

We now state and prove the main theorems of this section. We first show that the conditional hit- 
ting time of the continuous Euler process X given the discrete path observations Y tQ , . . . , Y tm may 
be written as an expectation of an explicit functional of the discrete path X tQ , . . . , X tm of the signal 
given the observations Y to , . . . , Y tm . Therefore, we reduce the initial problem to the characterization 
of the conditional distribution of X to , . . . , X tm given Y to , . . . , Y tm . 



Theorem 2.4. Set X m = {X to ,..., X tm ) and Y m = (Y to , . . . , Y tm ). We have: 

Ptf > tn\ Y tQ ,. . .,Y tm ) = *(Y to , . . . ,Y t J, 
where for y m := (y , . . . , y m ), 

m( y m ) = E[<S>(X m ) \Y m = y m ] . 
The function is defined for every z m = (zq, . . . , z m ) by 



m—l 



<s>(z m ) = F(t m ,t n , Zm ) rr c^y 1 ^ 

k 



with 



F(t m ,t n ,x) =E 



n-l 

riff 

,k=m 



k=0 



Xt k' Xt k+i 



(a)|X 



Proof. We may show similarly to Equation © in Lemma [2TT1 that 



P{t x >t n \Y t0 ,...,Y tm ) =E 1 { x>tm| F(t m ,t n ,^ m )|y t0 ,...,^ 



with 



F(t m ,t n ,X tm ) = P ( inf X u > a|X t „ 

V' il G[tm,tn] 



(15) 



(16) 
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It follows from Lemma [2731 that 



F(fmi 



E 



min X u > a|pQ fc ) fc=0 ,..., m ) \X tr 

u£[tm,tn] 



E 



n-l 



k=m 



n< fe ;f* +i (a)i^ 



On the other hand, combining a successive conditioning rule with the independence between W and 
W gives : 



P(r* >t n \Y tQt ...,Y tm ) =E E l {r x >tm} F(i m ,i n ,X tm )|(W tfe ,W tfe ) fe=0 



E 
E 
E 



E 



F(i m ,i n ,X t jP (rf > t m \(W tk ) k=0 



Xt m W[T* >t m \(x 



H k )k=0,...,m 



I Y m 

I Y m 
Y m 



Now, using Lemma [2731 once again, we have 

k=0,...,m 

so that 



m— 1 



fc=0 



fc CT fe 



(a) 



P(r a x > t n | y t0 , . . . ,Y tm ) = E [<D(X m ) I Y t0 ,. ..,Y tr 
where for every x m := (xo, . . . , x m ), 

m—l 

$ (x m ) = F(t m ,t n ,x m ) J] G^ +1 (a). 



fc=0 



This completes the proof. 



□ 



It remains now to characterize the law of x m \Y m = y m appearing in (fTBT ). This is the purpose 
of the following theorem, in which we shall write the conditional survival probability in a usual form 
with respect to standard nonlinear filtering problems. 

We set from now on a k = a(x k ,t k ),b k = b(x k ,t k ),5 k = S(y k ,t k ),h k = h(y k ,x k ,t k ), and v k = 
u(y k , t k ). We next give the main result of the paper. 



Theorem 2.5. We have: 

P(rf > t n \ Y t0 ,.. .,Y tm ) = $(Y to , . . . ,Y t J, 
where for y m := (y , . . . , y rn ), 

V[F(t m ,t n ,Xt m )K m L™] 



E[L 



(17) 



(18) 



with 



m—l 



y J 

m—l 



K m = ] [ G 

k=0 

The function g k is defined by 
9k(xk,yk',x k+ i,y k+ i 



a) and Ll 



I 9k(X tk 



(2irA k )Wj k 



exp 



k=0 



v k fx k+ i-m k y k+1 



m 



?x2 



25lA k V a k 



(19) 



with ml := x k + b k A k and ml := y k + h k A k . 
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Proof. It follows from Theorem [2j4] that 

P(rf > t n \ Y t0 , . . . , Y tm ) = E [<S>(X m ) \Y to ,..., Y tm ] 
where for every x m := (xq, . . . ,x m ), 

m—l 

<£ (x m ) = F(t m ,t n ,x m ) J] G^ +1 (a). 

k=0 

It remains to characterize the conditional distribution of X m given Y m . Recall that the dynamics of 
the processes (X tk ) and (Y tk ) are given for k = 0, . . . , m by 

[X tk+1 = X tk + b(X tk ,t k )A k + <7(X tfc ,t fc )(W ffc+1 " WtJ _ _ 

H + i = + /»(r th ,X tfc ,t fc )A fc + i/(lt fc ,tfc)(Wi fc+1 - W tk ) + 5(Y tk ,t k )(W tk+l - W tk ) 

(20) 

Let ~P XkL y k (X tk+1 G dx k+1 ,Y tk+1 G denote the density function of the couple (X tk±1 , Y tk+1 ) 

given (X tk ,Y tk ) = (x k ,y k ). Using Bayes formula and the Markov property of the process (X tk ,Y tk ) k 
we get: 

^xo,y {X tl G dxi, . . .,X tm G dx m \Y tl G cZyi, . . . ,Y tm G dy m ) 

_ p x , yo {{X tl G dxi, y tl G dyi); . . . ; (X tm G cfc m , V fm G cfy m )) 

_ F z ,yo(^i £ g rfyi) X • • • X Pa m -i,y m -i(^ m g <fam,*t m £ <fym) 

p xo,2/o( y *i £dyi,..., Y tm G dy m ) 

For every x m = (x , . - - , x m ) and y m = (y , y m ), set : 

N k (x m ;y m ) := V Xom (X tl G d Xl ,Y h G c%) x • • • x P :Em _ 1 , ?/m _ 1 (X tm G cte m ,y fm G dy m ). 
With this notation, the denominator reads: 

P 2 ,yo( y ti 6 • • • .^m G d 2/m) = J N k (x m ;y m )dx 1 . ..dx m . 

On the other hand, the random vector (X tk+1 , Y tk+1 ) | (X tk , Y tk ) = (x k , y k ) has a Gaussian distribu- 
tion with mean m k and covariance matrix T, k given for every fc = 0, . . . , m — lby 



m k 



Xk + t kA , k ), ^k = A k ( a * T\ 2 Y (2D 
y k + h k A k ) \v k v k v£ + °k) 



Then, for every fc = 0, . . . , m - 1, the density f k (x k , y k ; , y k+1 ) of (X tfc+1 , Y tk+1 ) | (X tfe , Y tk ] 
{x k ,y k ) reads 

tt \ 1 ( u k + 6 l ( (gfe±i - m i) 2 i (yfc+i - m lf 

f k {x k ,y k ;x k+1 ,y k+ i) = —— exp - 2 \ % 1 TT^2 

2u k 



(x fe+ i - m\)(y k+ x - m^)}^, (22) 



where m| := x k + 6fcAfc and m| := y k + h k A k . As a consequence, 

P (jf > t n \ (Y t0 ,. . . , Y tm ) = y m ) = E [<D(X m ) | (y, , . . . , Y tm ) = y m ] 

_ I <5>{x m )N k {x m ;y m )dx l ...dx m 
f N k {x m ;y m )dx 1 ...dx m 
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Now, we know that the random variable X tk+1 \X tk = x k has a Gaussian distribution with mean m\ 
and variance o\lS. k . Its density P k (x k , x k+ i)dx k+ i is therefore given by : 

Pk{x k ,x k+1 )dx k+1 = — exp( 1 2a fc jdxk+i- 

yJ2ir/\ k CF k V la k l\ k J 

Then, by definition of N k , we may write: 

N k (x m ;y m )dx 1 ...dx m 

m—l 



/ Y\ fk{x k ,y k ;x k+ i,y k+ i)dx k+ i 

J k=0 

[tj fk{x k ,y k ]x k+1 ,y k+1 ) 

f f'TJ fk(Xk,yk-,X k+1 ,y k+1 )\ T-T D , v , 

j f II 9k(x k ,yk 4 ,x k+ i,y k+ i) J P^^X^ G dxi,...,X tm £ cte m ) = E[L™], 



,fc=o / 

where is defined by (fl9l) and the next-to-last equality follows from the Markov property of the 
process (X tk , k = 0, . . . , m). Now, looking at the numerator of (|23T >. similar computations lead to: 



^{x m )N k {x m -y m )dx 1 ...dx m 

/m—l 
F(t m ,t n ,x m ) Y[ G^l^i^fkixk^k^Xk+i^k+^dxk 



I 



k=0 
m—l 

x k ,x k+1 



F(t m ,t n , x m) | ^ A k a k 



,k=0 



'm—l 



I 9k(x k ,y k ;x k+1 ,y k+ i) F xom (X tl e dxi,...,X tm £ dx r 



,k=0 



=V[F(t m ,t n ,Xt m )K m L™] 
Finally, going back to (|23T ). we obtain 



P K >t n \(Y t0 ,...,Y tm )=y« 



E[F(t m ,t n ,Xt m )K m Ll 



-.in 



which is the announced result. □ 

Remark 2.1. Remark that more generally we have for every real valued bounded function /, 

IE \f(X }L rn ~\ 

E[f(X tm )\Y t0 =y ,...,Y tM = y m ] = yi . (24) 

In this case we can estimate the right hand side quantity of (l24l using recursive algorithms. In our 
setting, we can adapt these algorithms by noting that the expression (fT8T ) may be read in the similar 
form of (EU) as: 

/ v \ E\F(t m ,t n ,X t )Z>1 

F > t n \ (Y t0 ,. . .,Y t J = y m ) = 1 1 E[L m ] 
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where 

m— 1 

L' y m = J] g k (X tk ,y k ;X tk+1 ,y k+1 )Gl tk f fc+1 (a) 



fc=0 



and where the involved functions are defined in Theorem [ 

Considering the model given in Example [TJ one may apply the result of Theorem 12.51 using the 
continuous Euler process X of the signal. However since in this framework the solutions of the 
stochastic differential equations are explicit, we refrain from using the Euler scheme to avoid adding 
additional error. In the next result we deduce a similar representation of the conditional survival 
probability using the explicit solutions of the stochastic differential equations of the model ©. 

Corollary 2.6. Consider that the signal process X and the observation process Y evolve following 
the stochastic differential equations given in (fJJ; 



dX t = X t {ndt + o-dWt), X = x , 

dY t = Y t {rdt + vdW t + 5dW t ), Y = y . 



(25) 



Then, 



W,\F(t m ,t n ,X t )K m L™] 
^>{y m ) = L V — mJ ^± (26) 



E[L 
with 

m—l m—l 

K m =U G*^f fe+1 (a) and L™ = ]J 9k (X tk , y k ; X tk+1 , y k+1 

k=0 k=0 

The function g k is defined by 



9k(xk,yk;x k+ i,y k+ i) = , — exp 



v 



2 ( log x k+ i - m\ log y k+l - m| n 2 



(27rA fc )3/2 7rfc *\ 25 2 A k \ a 

(27) 

with 

7r k := a 2 5x 2 k+l y k+1 , ml : = log x k + (ji - ^cr 2 )A k and m\ : = log y k + (r - ^u 2 - ^<5 2 )A fc . 
Proof. It follows from Ito formula that for every s < t, 

(X t = X s exp ((/i - \a 2 )(t - s) + a(W t - W s )) 



Y t = Y s exp ( (r - \v 2 - \b 2 )(t - s) + v{W t - W s ) + 5{W t - W s ) ) . 



(28) 



Then, for every k = 0, . . . , m - 1 the random vector (X tk+1 ,Y tk+1 )\(X tk ,Y tk ) = (x k ,y k ) has a 
bivariate lognormal distribution with mean m k and covariance matrix £ given by 

( \ogx k + ^-\a 2 )A k \ y- Al (° 2 ™ 

mk ~ \\og y k + (r - \v 2 - \5 2 ) A J ' ^ " Afc ^ v 2 + 5 2 

Hence, its density reads (setting fi k = aSx k +iyk+i) 

, , s 1 ( v 2 + 5 2 ( (log x k+1 - ml) 2 



2;/ 



+ 



CT^ 2 + 5 2 ) 

(logy fe+ i - m|) 2 



i/ 2 + J 2 



(logx fc+ i - m\) (log y fc+ i - m 2 ) 
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where m\ = log xj. + (fi — \cr 2 ) A& and m 2 k := log + (r — |v 2 — ^<5 2 ) Furthermore, for every 
fe = 0, . . . , m — 1, the random variable X tk+1 \X tk = x^ has a lognormal distribution with mean m\ 
and variance a 2 so that its density distribution Pk(xk, Xfc+i) is given by 



There exist several methods to estimate the above representation of the conditional survival prob- 
ability. These methods involve, amount others, Monte Carlo simulations and optimal quantization 
methods. Owing to the numerical performance of the optimal quantization method due for example 
to its fast performability as soon as the optimal grids are obtained, we will use optimal quantization 
methods to estimate the conditional survival probability. 

The use of optimal quantization methods to estimate the filter supposes to have numerical access 
to optimal (or stationary) quantizers of the marginals of the signal process. One may use the opti- 
mal vector quantization method (as done in the seminal work ifTTTl and used in 0) to estimate these 
marginals. This method requires the use of some algorithms, like stochastic algorithms or Lloyd's 
algorithm, to obtain numerically the optimal (or stationary) quantizers of the marginals of the sig- 
nal process. Given these optimal quantizers, the filter estimation is obtained quite instantaneously. 
However, the step of search of stationary quantizers is very time consuming. 

We propose here an alternative method, the (quadratic) marginal functional quantization method 
(introduce in lTT4ll to price barrier options), to quantize the marginals of the signal process. The 
marginal functional quantization method consists first in considering the ordinary differential equa- 
tion (ODE) resulting to the substitution of the Brownian motion appearing in the dynamics of the 
signal process by a quadratic quantization of the Brownian motion. Then, by constructing some 
"good" marginal quantization of the signal process based on the solution of the previous ODE's, we 
will show how to estimate the nonlinear filter. Since this procedure is based on the quantization of the 
Brownian motion and skips the use of algorithms to perform the stationary quantizers, the computa- 
tion of the marginal quantizers is quite instantaneous. This reduce drastically the time computation of 
the procedure with respect to the vector quantization method since we skip the step of the use of al- 
gorithms search of marginal quantizers by using instead, the marginals of the functional quantization 
of the signal process. 

In the rest of the paper we deal with the estimations methods of the conditional survival probabil- 
ity. 

3 Numerical tools 

Our aim in this section is to derive a way to compute numerically the conditional survival probability 
using Equation (PT8l) . To this end, we have to estimate three quantities: the quantity F(t m ,t n , X tm ), 
the expectation E[F(t m , t n , X trn )K m L™\ as soon as F(t m , t n , X tm ) is estimated, and finally, the ex- 
pectation E[L™]. Both expectations will be estimated using marginal functional quantization method. 
The probability F(t m , t n , X tm ) will be estimated by Monte Carlo simulations. 

Before dealing with the estimation tools, we shall recall first some basic results about both optimal 
vector quantization and functional quantization. Then, we will show how to construct the marginal 
functional quantization process which will be used to estimate the quantities of interest. 

3.1 Overview on optimal quantization methods 

The optimal vector quantization on a grid V = {x±, • • • , xn} (which will be called a quantizer) of 
an Revalued random vector X defined on a probability space (fi, A, P) with finite r-th moment and 




We then conclude the proof by using the same arguments than those of the proof of Theorem [23] □ 
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probability distribution Px consists in finding the best approximation of X by a Borel function of X 
taking at most N values. This turns out to find the solution of the following minimization problem: 

e N)r (X) = ml{\\X - X r || r ,T C R d , |T| < N}, (29) 

where \\X\\ r := (E[|X| r ]) 1 ^ r and where X r := ^^{XeCfr)} i s trie quantization of X (we 

will write X N instead of X r ) on the grid T and {Ci(T))i = i n corresponds to a Voronoi tessellation 
of R d (with respect to a norm | • | on F, d ), that is, a Borel partition of R d satisfying for every i, 

CAT) C \x G M d : \x - xA = min 

j=l,...,N J 

We know that for every N > 1, the infimum in ([29]) is reached at one grid T* at least, called a 
L r -optimal ^-quantizer. It is also known that if card(supp(Px)) > N then |T| = N (see e.g. [8] 
or [I^). Moreover, the L r -mean quantization error eN, r (X) decreases to zero at an ./V _1 / d -rate as 
the size N of the grid T goes to infinity. This convergence rate has been investigated in [CO and |[T5l 
for absolutely continuous probability measures under the quadratic norm on W 1 , and studied in great 
details in [H under an arbitrary norm on M d for absolutely continuous measures and some singular 
measures. 

From the numerical integration viewpoint, finding an optimal quantization grid V* may be a 
challenging task. In practice (we will only consider the quadratic case, i.e. when r = 2) we are 
sometimes led to find some "good" quantizations Proj r (X) (with T = {x±, . . . ,xn}) which are 
close to X in distribution, so that for every Borel function F : M. d \- > R, we can approximate 

E[F(X)] by 

N 

E[F(X N )]=Y,F(x i )p i , (30) 
i=i 

where pi = P(Proj r (X) = Xj). Amount "good" quantizations of X we have stationary quantizers. 
A grid T inducing the quantization X N of X is said stationary if 

V % £ j, xj and F (X £ U^dC^r)) = (31) 

and 

E \x\X N ] = X N . 

The stationary quantizers search is based on zero search recursive procedures like Newton algorithm 
in the one dimensional framework and some algorithms as Lloyd's I algorithms (see e.g. [0), the 
Competitive Learning Vector Quantization (CLVQ) algorithm (see |J5]) or stochastic algorithms (see 
lfT2l ) in the multidimensional framework. Note that optimal quantizers estimates of the multivariate 
Gaussian random vector are available in the website www. quantize . math-f i . com. 

We next recall some error bounds induced from approximating E[F(X)] by (1301) . Let T be a 
stationary quantizer and F be a Borel function on F, d . 

(i) If F is convex then 

E[F(X N )] <E[F(X)]. 

(ii) Lipschitz functions: 

- If F is Lipschitz continuous then (this error bound doesn't require the quantizer V to be 
stationary) 

\E[F(X)]-M[F(X N ))\ < [F] Up \\X-X N \\ 2 , 

where 

[F]Llp:=sup m^M. 

xj^y |-E y\ 
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- Let 6 : R d — > R + be a nonnegative convex function such that 9{X) G L 2 (F). If F is lo- 
cally Lipschitz with at most ^-growth, i.e. \F(x) — F(y)\ < [F]np\x — y\ (9(x) + 9{y)) 
then F(X) G L X (P) and 

|E[F(X)] - E[F(X N )]\ < 2[F] Up \\X - X N \\ 2 \\e(X)\\ 2 . 

(iii) Differentiable functional : 

If F is differentiable on M. d with an a-H6lder differential DF (a G [0, 1]), then 

\E[F{X)} - E[F(X N )]\ < [DF] a \\X - X N \\\ +a . 

Other error bounds related to the regularity of F may be found in |[T3"T0 . 

The optimal vector quantization may be extended to random vectors with values in a set of in- 
finite dimension, in particular to stochastic processes viewed as random variables with values in 
L 2 ([0, T],dt) endowed with the norm 



E ' 



2 



\X\ t2 



I E[X*]ds < +oo. 
Jo 

The functional quantization of the stochastic process (^t)t e [o,T] with dynamics 

dX t = b(X t ,t)dt + a(X t , t)dW t , X Q = xeR, 

is based on the functional quantization of the Brownian motion W. One way to quantize the Brownian 
motion is to use the optimal product quantization using its Karhunen-Loeve expansion which reads : 

n>l 



where £ n := (W,e n ) /\^K t , n > 1, is a sequence of i.i.d. random variables with standard normal 
distribution and 



/ \ 2 

1 



£n(t) = \ — sin 7r n — , A n = — : tit , n > 1. 

2 J T J ) ~ 

In fact, from the previous expansion, a functional quantization W of the process W of size at most N 
is defined by 

^ = E^fenW (32) 

n>l 

where |^ (with x Nn = {x± n , . . . , x 1 ^" }) is the optimal N n -quantization of £ n and N\ x • • • x N n < 
N, with Ni > 2 for i < n, and Nk = 1 for k > n + 1, so that the expansion defined in (l32l) is a finite 
sum. The product quantizer \ that produces the above Voronoi quantization W is defined by 

Xi(t) = ^2VK^Z n e n (t), i = (ii,--- ,w--) € II^ 1 ' ' iV ™}- 

ri>l n>l 

and for every multi-index 2 G rin>i {!> " ' ) -^n}> the associated Voronoi cell of x is 



n>l 
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The optimal product quantizer of size at most N, denoted W , of the Brownian motion is defined as 
the solution of the following optimization problem: 

mm{\\W - W N \\ 2 , iVi,--- ,N n > 1, Ni x • • • x N n < N, N> l}. (33) 

Moreover this optimal product quantizer induces a rate-optimal sequence of quantizers (see e.g. (9[ 
for more details), i.e. 

\\W - W N \\ 2 < K w 

(log N) 2 

for some real constant Kyv > 0. To define a functional quantization of the stochastic process (X t )t>o 
consider (x N )n>i a sequence of rate-optimal product quantizers of the Brownian motion and, for 
every multi-index i £ nn>i {!> " ' " > ^n}» with -/V*i X • • • x N n < N, consider Xj_ the solution of the 
following integral equation 

dxdt) = (b(xi{t),t) - ^aa'(xi(t),tf) dt + a(x i (t),t)d Xi (t), (34) 

where a' is the derivative of a. We define the functional (non-Voronoi) quantization X process of the 
stochastic process X of size at most N by 

N 

= j>a Ci(x) on n>i. 

i=l 

We have the following result. 

Proposition 3.1 (See fl9|). Under some suitable conditions on the coefficients of the diffusion, which 
in the homogeneous case are equivalent to: b is differentiable, a is positive twice differentiable and 
b' — b- — \o~o~" is bounded, we have 

\\X-X N \\ 2 = C((logiV)-i). 

We observe that for any initial value x of the quantized process X, the marginals X t are of size 
djv = Ni x • • • x N n where (N±, . . . , N n ) is the solution of (l33l . To define the marginal func- 
tional quantization process (still be denoted by (Xt)) of the process (Xt) we order the values of the 
marginals and define the marginal quantizer x dN (t) for every t > by x dN (t) = {x dN (t) , . . . , x d N (£)} 

andx djv (0) = x and define likewise the marginal quantizers of Wt by x dN (t) = {xf N (t)i ■ ■ ■ >X^(*)} 
The marginal functional quantization process of size at most N is then defined by 

d N 

X t = ^2 x i N ( t ) 1 {w t ec z (x dN (t))}- 

i=l 

Remai - k that the use of the marginal functional quantization method can not be justified from the 
theoretical point of view since we do not know yet the rate of convergence of the marginals of the 
quantized process to the marginals of the initial process. However, this method has proved its effi- 
ciency from the numerical viewpoint when used to estimate barrier option by optimal quantization, 
see |[T4l (not that the considered marginal quantization in |[T4l is a little bit different from the one con- 
sidered in this paper, nevertheless the numerical results are the same, up to at least a 1CP 3 absolute 
error order). 

Let us come back to the problem of interest and let the functionals Tv y m and w y m be defined for 
every bounded measurable function / by 

7r y , m f = V[f(X tm )L™K m } and w v , m f = E[f(X tm )L™]. 
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Then, 



Ily,mF(t m ,t n ,-) := P(r a > t n \ Y to = y , . . . ,Y tm = y m ) 



^y,m\ 



where l(x) = 1, for every real x. Then, it will be enough to show how to estimate 7r yim F(t m , t n , •) 
since w y:Tn l is estimated similarly. We discuss in the section below the estimation of ny im . 

3.2 Estimation of n y)Tn by marginal functional quantization 

Our aim is to estimate 

m—l 

TT y , m F(t m ,t n ,-) = E F(t m ,t n ,X tm ) Y[ 9k(X tk ,yk;X tk+1 ,y k+1 ) 

k=0 

by marginal functional quantization (with X = X in the model ©), where 

— — — — JC Xf 

9k( x t k , Vk]X tk+1 , y k+1 ) = g k (X tk , y k ; X tk+1 , y k+1 ) x * fc+1 (a). 

We deal with a general setting: the estimation of %,m/ for any bounded and measurable function 
/. Our main reference is ifTll where 7r yj7n f has been estimated using marginal vector quantization 
methods. Let us define for every k = 1, . . . , m the transition kernel H Vjk by 



#„,»/(*) = E[f(X tl )gZ_ 1 (x,y k - l ;X, t ,y t )\X tl _, 



X 



where P k (x, dz) is the density of the random variable X tk \Xt k _ 1 = x. We set 

H vja f = E[f(X Q )} = / f(z)fi(dz). 



(35) 



Then, we have (setting for every k, F x = a({X ti , i = 0, . . . , k})) 



k-l 



TTy :k f = E 



E 



E 



E 



E 



E 



f(X tk ) J] gf(X u , y t ; X u+1 , y i+1 ) | 

f{.X tk )9k-\{.Xt k _ 1 , Vk-i; X tk ,y k )\F t X k _^ Yl 9i(Xt i ,yi;X t . +1 ,y i+1 ) 

k-2 

f( x t k )g k _ 1 (X tk _ 1 , y k _x; X tk , y k )\X tk J gf(X u , yi]X u+1 , y i+1 ) 



k-2 



i=0 



where the last equality is a consequence of the Markov property of the process (X tk ). Thus, we 
deduce that for every k = 1, . . . , m, 



k-2 



Ky,kf = E Hy^ k f{X tk _ 1 ) Y[9i(X ti ,yi;X ti+1 ,y i+ i) 

= ^y,k-\Fly,kf ■ 

It follows that TTy^mf can be computed by the following recursive formula: 

Ky,mf = {Hyfl O FL y \ o • • • o Hy jTn )f. 



(36) 
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Therefore, to achieve the estimation of ix y ^ m , it remains to estimate the kernels H y ^- This will be 
done by marginal functional quantization. 

Consider time discretization steps t k , k = 0, . . . , m and let x Nk '■= {xi-> ■ ■ ■ > Xh k } ^ e a Nk- 
quantizer of Wt k (we will consider the marginal functional quantization of Brownian motion so that 
Nk = cIn for every k). Suppose that we have also access to the marginal functional quantization 
process (X tk )k of the process (X t )t>o over the time steps t k , k = m, . . . , n: {x\, . . . , x^ k }, of sizes 
N k (keep in mind that 7V& = for every k). 

It follows that the transition kernels H y k may be estimated for every k = 1, . . . , m by 

i=i 

where 

S^fe = 9k-i( x l-i^yk-r,^yk)Pk^ 1 = !>■■■ ,N k _r,j = l,...,N k (37) 
and where the p^'s correspond to the estimation of the transition probabilities from X tk _ 1 = x k _ l to 

K = 4- 

V l l = F(X tk = 4\Xt k -t = 4-l)> t = l,...,JV fc _i;j = l,...,JV fc . (38) 



Finally, one will perform the estimation TT y>m = H y fi o iT^ o • • • o -ff^m of n yjm as soon as we will 
be able to compute p\ J . In the proposition below we show how to compute these probabilities from 
the cumulative distribution function of the random variable W't fc |Wt fc _ 1 = x, which has a Gaussian 
distribution with mean x and variance (we refer to lTT4ll for a similar result). 

Proposition 3.2. The transition probabilities can be estimated for every k = 1, . . . , m by 

tfl « " N (x j k ~ ; xi-i) , (39) 

where jV(- ; xt_i) W the cumulative distribution function of the normal distribution with mean x\-i 
and variance and where for every k = 1, . . . , m, 

~xi + = ^;xr = ^ j = h-,N k 

1- N+ 

{ x k = -°°; x k = +°°- 

Proof. One has for every fc = 1, . . . , m, 

On the other hand, we have for every 

i , v ^ P(W tk < x;W tk , e Ciix^ 1 )) 

p(w tk < x\Wu , g cAx k - 1 )) = v k r fc ' 1 ; * — ,J -. 

Set Ci = Ci(x Nk ~ 1 )- Then the numerator on the right hand side of the previous equation may be 
expressed as 

P{W tk < x-W tk _ x G Ci)=J^ ( jf P(W ifc G dxlWi,., = i/JPCWi,., G <fr,))dx 

= / AA(x;y)P(Wt fc _ 1 G dy) (40) 
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where N(-;y) is the cumulative distribution function of Wt k \Wt k _ 1 = y. The last quantity is the 
approximation of (l40l > by optimal quantization with one grid point, considering that {xl-il is the 
quantizer of size one of the random variable Wt k _ 1 over the cell Ciix^- 1 ). □ 

Following the previous approach, the estimations Tf y , m and m y . m of Tr^m and w y ^ m are computed 
from the following recursive formulae: 



^,0 - Hyfi 



, m 



where 



and 



where 



H yk = fl£_i04-i, yk-r,xi,yk)Pk, i = l,...,d N ;j = l,...,d N ; 



w y$ — Yy,0 

^y,k ^y,k—l~^y,k • 



2^ii=\ - L y,fc a; y,fc-l 



,k = 1, . . . , m 



T vk = 9k-\(Xk-],yk-i;x{,y k )p^ , i = l,.. .,iV fe _i;j = 1, . . .,N k . 



Then we approximate n yjTn by IIj, )Tn given by 



A 1 ,, 



n^m — m b x 



i=i 



with 



n 



7T 



y,m 



y,m „JV, 



Z^j=l M y,k-l 



i = l,...,iV„ 



Recall that our aim is to estimate 



T^y,mF(trni tr 



where the function F{t m , t n ,-) is defined for every x > by 

"n-l 



n<:rf tfc+i (a)i^ m =x 



F(t m ,t n ,x) =E 
We shall take the estimation: 



(41) 



(42) 



(43) 



Remark that the initial function F(s, t, X s ) = P inf X u > a\X s (which has been estimated by 

\s<u<t J 

F(s, t, X s )) has semi-closed expression in some specific models, like in the model Q in which case 
it is given by 



inf X u > a|X s = NfaiX^t- s)) - — 

, s<u<t I \A. 



<7- 2 ( M -cr 2 /2) 



N(h2(X„t- 8 )) (44) 
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where 



hi(x, u) 



1 / 1 

—= y log(ac/a) + (fi- -a 2 )u 



h 2 {x,u) = — ■=( log(a/x) + U- \(? 2 )u) 



and where N(-) is the cumulative distribution function of the standard Gaussian distribution. 

Except in these specific cases, the function F(t m ,t n , •) has to be estimated, and we shall do it by 
Monte Carlo methods. 

3.3 Estimation of F(t m , t n , ) by Monte Carlo 

It remains to estimate the function F(t m: t n , x) defined by 



F(t m ,t n ,x) = E 



n-l 

no 

,k=m 



X tk' Xt k+t 



fc°fc 



(B)\X im =X 



by Monte Carlo simulations. The steps of the Monte Carlo procedure are the following. 

1. Let us consider regular time discretization steps {t m , £ m +i, . . . , t n } over [t m , t n ] and let M be 
the number of trials. We simulate for every j = 1, . . . , M the sample path (Xf . )fc=m,...,n» with 
X\ = x for every j. 



2. Setting 



we estimate F(t m ,t n ,x) by 



n-l 



k=m 



a) := n G £f fc+1 ( a 



M 



F (trmt n ,x) — ^ 'y ] fP m n (x, a). 



(45) 



Consequently, integrating both formulae (1431 and (l45l) . the conditional survival probability 

will be estimated (for a fixed trajectory (yo, • • • , y m ) of the observation process (Y to , . . . , Yt m )) by 



(46) 



t=l j=l 

Remark 3.1. It follows from Monte Carlo error analysis that for every x > 0, 



1 



<M 



(47) 
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3.4 Error analysis 

We target to give in this section the error bound resulting from the estimation of 

F(T*>t\Y t0 ,...,Y tm ) 

by 

N m M 

n^m-F (t m ,t n , •) = — y ]y ^ ^y,mPm,n( x m'i a )- 

i=l j=l 

We shall first give the error bound due to the fact that we compute the survival function of the random 
variable t x instead of t x , and then the error bound due to the simulations procedures. To this end, 
the following definitions and assumptions are needed. 

Definition 3.1. A probability transition P on E is C-Lipschitz (with C > 0) if for any Lipschitz 
function / on E with ratio [f]up, Pf is Lipschitz with ratio [Pf]up < C[f]up- Then, one may 
define the Lipschitz ratio [P]it P by 

[P]up = sup { ^\J} Lip , f a nonzero Lipschitz function } < +oo. 
\t\Lip 



Remark that that in our framework, the transition operators Pk(x, dy), k = 1, ■ ■ ■ ,m are Lips- 
chitz, so that we set 

[P]Lip ■= max [P k ]Li P < +oo. 

fc=l,"- ,m 

We furthermore define some useful quantities which appear in the error bound in the following. 

(i) For every k = 1, • • • , m, we set 

K" 1 := max ||g fe ||oo, 
a fc=l,— ,m 

where Halloo is the supremum norm of the functions g% defined by (fT9l ). 

(ii) For every k = 1, ••• ,m, let [gl]up and [g^Lip be so that for every x,x',x,x' G E and 

bl(^) y 5 z'> y') -fK^y^jyOI ^ [fi'fclLiply,^) \x-x\ + [gl]Li P (y,y') W - x'\. 

Let us make now some assumptions which will be used to compute (see Q) the convergence rate 
of the quantity E|lr x >t \ — l{ T >t}| towards 0. 

(HI) 6 is a C£°(R, R) function and a is in C 6 °°(R, R). 

(H2) there exists do > such that Vx G R, a(x) 2 > <Tq (uniform ellipticity). 

Before giving the error bound associated to our estimation we recall the following useful results. 
Consider in this scope that 

t x = \ni{u > 0, X u D} 

where D = (a, +oo) and X is the signal process. Let (X tk )k=o,...,m tne continuous Euler process 
taken at discrete times t^, k = 0, . . . , m and 

t x = inf{n > 0,X U D}. 

We have the following result. 
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Proposition 3.3 (see Q). Let t > 0. Suppose that Assumptions (HI) and (H2) are fulfilled. Then, 
for every t] £ (0, i[ there exists an increasing function K(T) such that for every t G [0, T] and for 
every x G E, 

1 if (T) 

where n is the number of time discretization steps over [0, t}. 



l 1 ^^} - 1 {r^>t}l 



< 



The convergence rate of the filter approximation is given by the following theorem. Since we do 
not know the convergence rate and some properties as the stationary property of the marginals of the 



functional quantization we consider here that for every k, X tk k denotes the marginal quantization of 



(X tk ) of size N k , obtained from the marginal vector quantization method (see ifTTTO . 
Theorem 3.4. (see HI IV for a similar result). We have for every p > 1, 



\U y , rn F(s, t, •) - U y , m F(s, t,-)\< 9 ]T Bf(F{s, t, -),y,p) \\X tk - \\ p 

<f>m{y) V0 m (y) k=0 



where 

<t>m{y) ■= ^y,ml, 4>m(y) ■= 7Cy,ml 

and where 



B?(/,y,p) := (2 - S 2 , p )[P]^[f] Lip + 2 ^j^([gl +1 ] Lip (y k ,y k+1 ) + [g 2 k ] Lip (y k ^,y k )) 
+ (2-S 2 J^ J2 [P]i^" 1 (b)] Ljp ( yj _ 1 , % + [ J P]L, P bJ]L, P (%-i,y J ))) 

9 j=k+l ' 
(with the convention that g§ = g m +i = and 6 niP is the usual Kronecker symbol). 
Proof. The proof is similar to the proof of Theorem 3.1 in iTTTTl . 

□ 

Remark 3.2. One may remark that the function B™(/, y,p) involves the norm ||/||oo which for 
/ = F(s, t, •) is bounded by 1. 



Let us give now the error bounds induced from the approximation of P(t* > t n \Y to , . . . ,Y t 



by Tly tm F (t m , t n , •). 



Theorem 3.5. (See /El/ far a similar result) Suppose that the coefficients b and a of the continuous 
signal process X are such that Assumptions (HI) and (H2) are satisfied and let r] G (0, i]. Then, 



1 N m M 

^{ T a. > tn\Yt , ■ ■ ■ ,Yt m ) — — 'y ^ li vr ^ mn (x m , a 



N m M 



M 

i=l j=l 



+ - : r-S : : £ B ™(^ *> 0, y. p) " *f * Hp, 
9m{y) v ^m(y) fc=0 



where the functions (f) m , <f> m , B™(-, are introduced in Theorem \3.4 
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Proof. We have 



N m M 



P^ 7 "^ > M^toJ • • • j ^m) ~~ ^Kj,mPrn,n( x mi a ) 

i=l j=l 

< |P(r* > t n |y i0 , . . . ,y t J - P(r* > t n \Y t0 , 

I rfj/jm-^X^m , t n , •) Ily ^rnF{tm i^ni ') | 



> ) | 



(48) 



Owing to RemarkOand to the fact that H* m < 1 we get 



iv„ 



\^-y,mF(^S,t, n.y. m i^ (tmitni 



i=l 



i=l 



iV„ 



m ,t n , X) 



i=l 



x>0 



< sup |F(t m ,t n ,x) - F M (t m ,t n ,x)| 

a;>0 

On the other hand, the error bound of the term \Uy tm F(t m ,t n , •) — Uy )m F(t m ,t n , -)| is given by 
Theorem 13.41 

Now, let us consider the first term of the right hand side of Equation (l48l) . We have: 

|lP(7a X > t n \Y t0 , ...,Y t J- P(rf > t n \Y t0 , ...,Y tm ] 



E 



= E 
< E 



E 



1 



11 



{*/>*»} ~ 1 {T*>t n } 



the last statement following from Proposition [ 



□ 



3.5 Numerical examples 

We deal with numerical simulations in this section by considering two example of models. 
In both examples we fix t m = 1 and, given a (simulated) trajectory of the observation process Y from 
to t m , we estimate the conditional cumulative function P(r^ > t n \Y to , . . . ,Y tm ) using formula 
(l46l) . for t n varying 0.1 by 0.1 from 1.1 to 11 (where the time unit is expressed in years). Furthermore, 
we set the number m of discretization points over [0, t m ] equal to 50 and for every k = 1, • • • , m, 
the quantization grid size is set to 966 (as a consequence of the numerical solution of the Problem 
[33] for N = 10000, with the optimal decomposition (23,7,3,2), see ifTBl for more detail), with 
Nq = 1. All the programs have been coded using the C language on a CPU 2.7 GHz and 4 Go 
memory computer. 

Example 3 (The "Black-Scholes" example). 

The first model is the one considered in Example [J and Corollary 12.61 where the dynamics of the 
signal process X and the observation process Y are given by 



dX t = X t (jidt + adW t ), X = x , 

dY t =Y t {fidt + adW t + 5dW t ) , Y = y 



(49) 
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or equivalently 



dYt 



— = -r^ + SdWt- 



We choose the following parameters (as in Q) of the model: \i = 0.03, a = 0.03, xq = yo = 
86.3 and a = 76. The numerical results are depicted in Figure Q] and Figure |2 In Figure [T] we 
draw three trajectories of the observation process Y for the same 5 = 0.1 (left side graphic) and 
the corresponding cumulative functions P(r^ < t n \Y to , . . . ,Y tm ) on the right hand side graphics, 
t m = 1 and t n G [1.1, 11], in years. We remark that, for a fixed time t n , the lower the trajectory is, 
the higher its probability is to hit the barrier. 

The left hand side graphics of Figure|2]corresponds to three trajectories of the observation process 
Y for 5 G {0.1, 0.3, 0.5} and the right hand side graphics, to (a zoom of) the corresponding cumula- 
tive functions P(r^ : < t n \Y to , . . . , Y tm ), with t m = 1 and t n G [1.1, 11], in years. We observe in this 
example that the noisier the observations are, the higher the probability is to hit the barrier a before a 
fixed time t n . 

Note that in both examples, the function F{t m ,t n: •) has been computed using formula (1441) and 
the computation time to get one cumulative function P(t^ < t n \Y to , . . . , Y tm ) for a given t n is about 
of 24 seconds. 



U51 
110^ 



Y UP: delta = 0.1 

Y MID: delta = 0.1 

Y DOWN: delta - 0.1 




0.1 0.2 0.3 



0.5 0.6 0.7 0.8 0.9 1.0 




Figure 1: ("Black Scholes example") Three trajectories for the observation process Y and for S — 0.1 (on the left), and 
the corresponding cumulative functions P(r^ f < t n \Yt ,Yt m ) with t m = 1 and t n G [1.1,11] years (on the right). 



Example 4 (The Ornstein-Uhlenbeck example). 

In the second model, we suppose that both the signal and the observation process evolve following 
the Ornstein-Uhlenbeck dynamics: 




- X t )dt + adWu X = x , 

Y t )dt + adW t + 5dW t , Y = y 



(50) 



or setting Z t = Y t - X t 



dZ t 



-XZ t dt + SdW t . 



meaning that Z is still an Ornstein-Uhlenbeck process with mean value 9 = and with volatility 5. 
The parameters are chosen as follows: A = 0.18, 6 = 0.35, a = 0.12, xq = yo = 0.35 (as in (3l) 
and a = 0.2. The numerical results are represented in Figure [3] where we depict three trajectories 
of the observation process Y for 5 = 0.16 (left hand side graphics of Figure [3]) and the associated 
cumulative functions P(t^ < t n \Y to ,. . . , Y tm ), with t m = 1 and t n G [1.1, 6] in years (right hand 
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Figure 2: ("Black Scholes example") Left graphics corresponds to three trajectories for the observation process Y, for 
S = 0.1 (Y UP), 8 = 0.3 (Y MID), S = 0.5 (Y DOWN), and the right hand side graphics correspond to a zoom of the 
cumulative functions P(t* < t n \Yt , . . . , Yt m ) withi m = 1 and t n G [1.1, 11] years. 



side graphics of Figure[3]). Once again, we remark that, as in the "Black Scholes example", for a fixed 
time t n , the lower the trajectory is, the higher its probability is to hit the barrier. 

In this example, the function F(t m ,t n , •) has been computed using Monte Carlo simulations of 
size M = 10 5 (see the formula (|45T>) with 50 discretization steps over [t m ,t n ]. The computation time 
to get one cumulative function P(t;^ < t n \Y to , . . . , Y tm ) for a given t n is about 6.5 minutes. 




Figure 3: ("Ornstein-Uhlenbeck example") Left graphics correspond to three trajectories for the observation process Y, 
for 5 = 0.16 and the right hand side graphics correspond to the cumulative functions P(r^ f < t n |Y to , . . . , Y tm ) with 
tm = 1 and t n G [1.1, 6] years. 
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